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We report on a search for the X(1812) state in the decay — > K^uKp with a data sample of 
657 X 10® BB pairs collected with the Belle detector at the KEKB e'^e collider. No significant 
signal is observed. An upper limit ii'±X(1812), X(1812) uj4>) < 3.2 x 10"^ (90% C.L.) 

is determined. We also constrain the three-body decay branching fraction to be B{B^ — > K^ujtj)) 
< 1.9 xlO^'' (90% C.L.). 

PACS numbers: 12.39.Mk, 13.20.He 



Using a sample of 5.8 x lO'' J/tp events, the BES collab- 
oration observed a near-threshold enhancement in the uj(f> 
invariant mass spectrum from the double OZI suppressed 
J/if) ^ 7W^decay with a statistical significance of more 
than lOcr [l|]. When fitted with a Breit-Wigner, this 
enhancement, called X(1812), has the following mass, 
width, and product of branching fractions: 

M = (18121^^ ± 18) MeV/c2, 
r = (105 ±20 ±28) MeV/c2, 
S( J/i/; -.^X,X^ oj(j)) = (2.60 ± 0.27 ± 0.65) x 10'^. 

Partial wave analysis favors a spin-parity assignment 
of .F^ = 0++ for the X(1812). In the related u^j 
mode, Belle has seen a dramatic threshold enhancement 
in B+ K+ujip, the Y(3940) 0], which has now been 
confirmed by BaBar 0]. 

If the X(1812) is a qq meson, the X(1812) uj(j) 
branching fraction should be very small due to OZI sup- 
pression and the limited available phase space, in con- 
trast with the BES observation. Suggestions have been 
made that the X(1812) may be a tetraquark state (with 
structure Q^Q^), since some tetraquark states decay to 
vector-vector mesons dominantly by "falling apart" and 
their masses are at the threshold of two vector mesons Q . 
Other works speculate that it may be a hybrid [1| , glue- 
ball state , an effect due to intermediate meson rescat- 
terings or a threshold cusp attracting a resonance Q . 
In this paper, we report our search for this state in the 
decay K^u)(f>. On the other hand, this decay pro- 

ceeds via a b ^ s penguin with ss and uu popping. A 
similar decay mode _B+ K^(j)<j), which proceeds via a 
b s penguin diagram with double ss popping, is the 
only observed charmless B VVP (two vector mesons 
and one pseudoscalar meson) mode and has a rather large 
branching fraction [(4.9t|^) x 10"^] [i,[l^. Therefore, 
even if the X(1812) cannot be observed, measurement 
of the B^ K^LOcj) three-body decay is also helpful for 
investigating decay mechanisms. 



This analysis uses 605 fb ^ of data containing 657 x 
10^ BB pairs. The data was collected with the Belle 
detector flil at the KEKB 



]J e+e- 



asymmetric-energy 
(3.5 GeV on 8.0 GeV) coUider operating at a center-of- 
mass (CM) energy of the T(4S') resonance. 

The Belle detector is a large-solid-angle spectrome- 
ter [ll[. It consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Chcrcnkov counters (ACC), time-of- 
flight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals located 
inside a superconducting solenoid that provides a 1.5 T 
magnetic field. An iron flux return located outside the 
coil is instrumented to detect K'^ mesons and to identify 
muons (KLM). 

-S-daughter candidates are reconstructed from the de- 
cays Lu 7r"'"7r~7r° , (f) — > K~^K~ and 7r° — ^ 77. Charged 
tracks are identified as pions or kaons by combining infor- 
mation from the CDC, ACC and TOF systems. We re- 
duce the number of poor quality tracks by requiring that 
\dr\ < 0.3 cm and \dz\ < 1.5 cm, where \dr\ and \dz\ are 
the distances of closest approach of a track to the interac- 
tion point in the transverse plane and z direction (oppo- 
site to the direction of the positron beam), respectively. 
In addition, tracks matched with clusters in the ECL 
that are consistent with an electron hypothesis are re- 
jected. We use a kaon identification likelihood ratio Rk.tt 
= Lk /{Lk+Lt^) to discriminate K and tt candidates. 
The requirements i?A',7r>0.4 for a kaon and i?^^ 7r<0.6 for 
a TT are used. The efficiency to identify a kaon(pion) is 
94%, while the probability that a pion(kaon) is misidenti- 
fied as a kaon(pion) is about 10%. Candidate 7r° mesons 
are reconstructed from pairs of photons, where the en- 
ergy of each photon in the laboratory frame is required 
to be greater than 50 MeV. We select 7r° mesons with an 
invariant mass in the range 0.1193 GeV/c^ < M(77) < 
0.1477 GeV/c^ and a momentum in the laboratory frame 
P^^*^ > 0.38 GeV/c. 
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Particles satisfying the above selection criteria are then 
used to reconstruct w and (j) mesons. We select can- 
didates in the invariant mass windows 0.75 GeV/c^ < 
M^+^-^o < 0.81 GeV/c2 and 1.00 GeV/c^ < Af^+x- < 
1.04 GeV/c^. A vertex fit for the </> and uj candidates is 
also performed. In addition, we require three kaons in 
the final state, one directly from the S-meson decay and 
the other two from the (j) decay. To distinguish the two 
kinds of kaons and reduce multiple candidates, we require 
kaons from the cf) to have momenta pk± < 1-5 GeV/c in 
the CM frame. 

Candidate K^uxf) decays arc identified by us- 

ing the energy difference (AS) and the beam-energy- 
constrained mass (Mbc). These are defined as AE = 

Eb - -Bboam and Mbc = y/^l^i^n^-PB^ whcrC E^beam 

denotes the beam energy, Eb and pb denote the re- 
constructed energy and momentum of the candidate B- 
meson, all evaluated in the e+e^ CM frame. We se- 
lect events satisfying \AE\ < 0.2 GeV and 5.20 GeV/c^ 

< Mbc < 5.29 GeV/c^, and define signal regions 
-0.15 GcV < A£; < 0.05 GeV and 5.27 GeV/c^ < Mbc 

< 5.29 GeV/c2. 

The dominant source of background arises from ran- 
dom combinations of particles in continuum e~^e^ qq 
events {q=u,d,s,c). To discriminate spherical BB events 
from jet-like qq events, we use event-shape variables: 
specifically, 16 modified Fox- Wolfram moments [l3| com- 
bined into a Fisher discriminant, T [l3|. Additional dis- 
crimination is provided by 9b, the polar angle in the CM 
frame between the B direction and z direction. Correctly 
reconstructed S-mesons follow a (1 — cos'^6'b) distribu- 
tion, while fake candidates from continuum tend to be 
uniform in cohOb- 

Further continuum background suppression is achieved 
using 6-fiavor tagging information. The Belle flavor tag- 
ging algorithm [15| yields the flavor of the tagged meson, 
^(=±1), and a flavor-tagging quality factor, r. The lat- 
ter ranges from zero for no flavor discrimination to one 
for unambiguous flavor assignment. For signal events, q 
is usually consistent with the flavor opposite to that of 
the signal 5, while it is random for continuum events. 
Thus, the quantity qrFs is used to separate signal and 
continuum events, where Fb is the charge of the signal 
B: Fb =+1 (-1) for B+ (S") 

We use a Monte Carlo (MC) sample [l^ to form T 
and to obtain the cosOb and qrFB distributions. Proba- 
bility density functions (PDFs) are derived from J- and 
the cosOb distributions and are multiplied to form signal 
(Lg) and continuum background {Lqq) likelihood func- 
tions, which are further combined to form a likelihood 
ratio Rs = Lg/ {Ls + Lqq). We divide events into six qrFB 
bins and determine the optimum selection criteria for 
each bin by maximizing Ng/^/Ng + Nb, where Ng is the 
number of signal MC events in the signal region, and Nb 
is the number of background events estimated to be in the 
signal region by assuming B{B^ K^ojcf)) = 1.0 x 10~^. 



This optimization preserves 57.9% of the signal while re- 
jecting 98.6% of the continuum background. 

Applying all of the above criteria, the fraction of events 
having multiple candidates is 21%. To select the B- 
mcson candidate, we add the of the lo meson vertex flt, 
and the of a tt" ^ 77 fit constrained to the PDG [l3| 
value of 7r° mass: the candidate with the smallest value 
is chosen. If multiple B candidates still remain, we use 
the of the (j) meson vertex fit and of the 5-meson 
vertex fit to choose the best one. 




FIG. 1: Projection of the data (points with error bars) and 
fit results onto (a) Mbc, (b) A£, (c) M^+^-^o, (d) Mj^+j^^ 
with the other variables satisfying Mbc £ (5.27, 5.29) GeV/c^, 
AE e (-0.15,0.05) GeV, M^+^-^o £ (-0.75,0.81) GeV/c^ 
Mk+k- gJI.OO, 1.04) GeV/c^ : signal (dot-dashed), qq 
(dashed), BB (dot-dot-dot-dashed) and total (sohd). 

In addition to the dominant continuum background, 
charmed B decay (6 c) and charmless B decay 
(6 u, d, s) backgrounds are studied using dedicated 
MC samples that are respectively about two and 25 
times the size of the data sample. Charmless B decay 
background is found to be small and is neglected. The 
following charmed B decay channels are studied using 
dedicated Monte Carlo samples: Dfn^TT^ 
Df ^jr±(/> ; B± ^ DJtt^K^, Df 7r°7r±(/) and 
B± ^ D^jD^K^, D^jD^ ^ Ti+TT-TT^K+K-. To mea- 
sure the three-body B^ K^ujcf) branching fraction, 
we require Mx+k-tt+tt-tt" > 2.2 GcV/c^ to exclude 1?° 
background; and require {M^^i^+x- — to_dJ>0.15 GeV/c^ 
as well as \M^Oy^x+K- — yran J>0.15 GeV/c^, where 
is the nominal Dg mass [17|, to veto the Dg background. 

We obtain the signal yield using a four-dimensional 
extended unbinned maximum likelihood (ML) fit to AE, 
Mhc, Mtttttt and Mkk- The likelihood function consists 
of the following components: signal decays, continuum 
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background {qq) , and charmed B-decay background (6 — > 
c). For all components, no sizable correlations are found 
among the fitting quantities. The PDF for event i and 
component j is defined as 

V] = V.iAE^) X V,{MU X r,{Mi,j,) X V,{Ml^^). (1) 

The signal Afbc is parameterized by the sum of a sin- 
gle Gaussian and an ARGUS function [l^, AE by a 
Gaussian and a Crystal Ball function [l^, and A/^+a'-; 
M^+^-ttO by Brcit-Wigner functions. For continuum 
background, Afbc is parameterized by an ARGUS func- 
tion, AE by a second-order Chebyshev polynomial, and 
Mk+k- : -^ir+ir-TrO by the sum of Breit-Wigner functions 
and first-order Chebyshev polynomials. BB background 
modelling is similar, but with a first-order Chebyshev for 
AE and for Mt^+^-.^o. All function parameters are de- 
termined from MC simulation. 

In our final fit to the data, the signal and qq yields are 
allowed to vary; the fraction of & — > c events is very small 
and thus the yield is fixed in the fit according to MC. 

The likelihood function to be maximized is given by 



£ = 



"Iv! 



N 



(2) 



where Yj is the yield of events for component j and N is 
the total number of events in the sample. 

Figure [T] shows the fit results. Peaking behavior ob- 
served in AE, Mbc, AfTTTTTT and Mkk is consistent with 
that from MC expectations. The branching fraction 
is evaluated using the following quantities: the signal 
yield Y^^k = 22.1^7 2 with reconstruction efficiency 
£ ~ 7.04 X 10~^: the combined daughter branching frac- 
tion Sd ~ 0.439 [l7|; a correction of 0.946 to the efficiency 
of A'/tt identification requirements, which takes into ac- 
count small differences between MC and data; and a total 
of 657 X 10® produced BB pairs, where equal fractions of 
B^B~ and B^B^ are assumed. 



The sources of systematic error are listed in Table HI 
The quoted 6% track reconstruction efficiency is from 
the consideration that there are five tracks in a selected 
event and for each track the efficiency error is 1.2%. The 
errors due to continuum suppression requirements are ob- 
tained by varying these cuts while the errors on the PDF 
shapes are obtained by varying all fixed parameters by 
±1(7. Toy MC tests and GEANT-based Detector Simula- 
tion (GSIM) tests are performed, we find that the fit bias 
can be neglected. To estimate the error due to the b ^ c 
contribution, we vary the normalizations by ±50%. 

Our final result for the three-body branching fraction 
based on the 605 fb~^ data sample is 



B{B^ K^Lo4>) = (1.15 



+0.43+0.14 
-0.38-0.13 



) X 10' 



where the first error quoted is statistical and the second 
systematic. We obtain the 90% confidence level upper 



limit B{B^ 



< 1.9x10 °bya frcquentist 



method using ensembles of pseudo-experiments. For a 
given signal yield, 10000 sets of signal and background 
events are generated according to the PDFs, and fits are 
performed. The confidence level is obtained from the 
fraction of samples that give a fit yield larger than that 
of data (22.1). We take into account systematic errors 
by varying the fit yield by the total systematic errors 
described in Table [D The significance of the signal, esti- 
mated using this method, is 2.8(T. 




2.2 2.3 



TABLE I: Systematic errors for ^ K'^LOfj)). In cases 

where the error on B(B± -> 7<:±X(1812), X(1812) ^ LU(j>) is 
different, it is shown separately in parentheses. 



Type 


Fractional 


error (%) 




+ (T 


— O" 


Tracking 


6.00 


6.00 


K/tv id 


2.90 


2.90 


it" Reconstruction 


4.00 


4.00 


Daughter B 


1.45 


1.45 


Signal/Background Modeling 


5.25(31.1) 


2.50(22.7) 


BB Background Yield 


1.75(2.10) 


0.89(1.60) 


MC Statistics 


0.56(1.06) 


0.56(1.06) 


Continuum Suppression 


7.21(12.0) 


7.26(12.4) 




1.36(1.36) 


1.36(1.36) 


Total 


12.1(34.4) 


11.2(27.1) 



FIG. 2: Mass spectrum in the ujcj} fit with the following 
components; B"*" — > Kuxf) three-body (dotted), BB (dot- 
dashed), qq (dashed), (dot-dot-dashed), _Ds (dot-dot-dot- 
dashed), B± i«:±X(1812)(long-dashed), and total(solid). 
The spectrum is also shown in the inset with an expanded 
vertical scale 

We next study the iu(f) mass spectrum. Because the 
aforementioned Mj^+j^-t^+t^-^o , lA/yrK+if- — "^nj and 
|Af^o^^+^- —rriD^l niass cuts influence the shape of the 
LU(j} invariant mass spectrum, we did not use them and fit 
the and Dg backgrounds simultaneously. 

We produced 0.6 million (2.0 million) (£»,,) 
background MC events for the decay B^ 
fs:±A(1812), A(1812) ujcf). The A(1812) mass and 
width are taken from the BES measurement and its 
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PDF is modeled by an rARGUS (reversed ARGUS [l8|, 
^'rARGUs(a;) = -FARGUs(2t - x), where t is the thresh- 
old) function plus a Breit-Wigner with a threshold. The 
three-body decay PDF is an rARGUS function, the 
background PDF is the sum of an rARGUS function and 
a Breit-Wigner, the Dg background is the sum of an rAR- 
GUS function and a Gaussian, while the qq, BB back- 
grounds are also modeled by rARGUS functions. We ob- 
tained all the parameters from MC samples. In our final 
fit to the data, we fixed the yield of and Dg back- 
grounds according to the PDG branching fractions 17[, 



and fixed the yield of BB background according to MC 
simulation. 

The final result is shown in Fig. [51 No significant 
signal is observed; the yield of the X(1812) is Q.2'^\ 
events. The systematic errors are also listed in Table. U 
where those in parentheses are for the items that differ 
from those in the three-body decay analysis. We also in- 
clude the errors from the fraction of , Dg background 
and the X(1812) width into signal/background modeling. 
Using the pseudo-experiment method described above 
and taking the systematic errors into account, wc find 
a limit on the product branching fraction of B{B^ 



is:±A:(1812),X(1812) 



< 3.2 X 10-^ (90% C.L.) 



In summary, using a data sample of 605 fb~ col- 
lected with the Belle detector, we present a search 
for the X(1812) meson in the decay B^ K^ujcj). 
No significant signal is observed. An upper limit for 
the product B{B^ if±A:(1812), X(1812) ^ luc/)) < 
3.2 X lO"'^ (90% C.L.) is determined. We also mea- 
sure the three-body B^ — s- K^ojcj) decay branching frac- 



tion 6(B± 



[1.15 



-1-0.43 -1-0.14 
-0.38 -0.13 



(< 1.9)] X 10' 



where the upper limit is at the 90% confidence level. 
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